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We study the effects of the addition of the charm quark on the QCD equation of state at zero and 
nonzero chemical potential on lattices with N t = 6. Our ensembles are quenched with respect 
to charm and the charm quark is a valence staggered quark. Along the trajectory of constant 
physics the ratio m s /m c is kept constant after tuning the charm quark mass at a lattice spacing of 
about 0.09 fm. We find that the charm quark has a significant contribution to the equation of state 
at zero chemical potential already at temperatures between about 1.2T C and 2T C . The additional 
contribution at nonzero chemical potential vanishes within the current statistical uncertainty. 
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1. Introduction 

The equation of state (EOS) of the quark-gluon plasma is an essential input to the hydrody- 
namics models which are used to interpret the experimental data from heavy-ion collisions. The 
experiments at RHIC create a fire ball which thermalizes within x « 10~ 24 s [jl]]. The u, d and s 
quarks participate in the thermal ensemble describing the state of the thermalized fire ball. On the 
other hand, under the experimental conditions the c quark probably is not thermalized and thus the 
2+1 flavor EOS is considered sufficient for the hydrodynamics models applied to the experimen- 
tal data. However, a quark-gluon plasma existed microseconds after the Big Bang. Under these 
primordial conditions the c quark probably participated in the thermal ensemble as well, which 
implies that for the study of the early Universe, the EOS with 2+1+1 flavors might be important. 
Previously, the question of the charmed quark contribution to the EOS at zero chemical potential 
has been studied in [Q]. That work was performed with the p4 action on lattices with N t = 4. 

In the present study, we determine the charm quark contribution to the EOS at zero and nonzero 
chemical potential. The calculations are done on the MILC collaboration 2 + 1 flavor asqtad lattices 
with mi = 0. lm s , where m s is tuned close to the strange quark mass. The nonzero temperature 
calculations are done on lattices with N t = 6. We simulate on trajectories of constant physics. 
Similarly to Ref. [Ep, the charm quark is just a valence staggered quark. The charm quark mass 
is tuned to match the rest mass of the D s at j8 = 7.08 {a « 0.086 fm) where the discretization 
effects are smallest on our trajectory. Within 4% at that point, m c /m s = 10. We have kept that ratio 
constant for lower temperatures. The tuning is probably incorrect at the lowest temperatures, but 
we do not expect this to matter much due to the large mass of the c quark and its small contribution 
there. The nonzero chemical potential calculation is done using the Taylor expansion method [^] 
to sixth order. We used 800 random sources in the transition region and 400 outside it. For the 
calculation at nonzero chemical potential, the valence c quark has a low cost in terms of computer 
time, but it requires a sizable software development. For 2 + 1 flavors we have 95 observables to 
code and for 2 + 1 + 1 flavors there are 399. 

We determine the EOS at zero chemical potential using the integral method [|], ||], in which 
the pressure is obtained from the interaction measure (trace anomaly) by integrating it along a 
trajectory of constant physics. At nonzero chemical potential, the Taylor expansion method gives 
for the pressure 



where fXi^ c are the chemical potentials in physical units for the light (u,d), heavy (s) and charm 
(c) quarks. Due to CP symmetry the terms in the above are nonzero only if n + m + k is even. The 
coefficients are 
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A similar expression can be written for the interaction measure coefficients. 
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Figure 1: Interaction measure (/), pressure (p) and energy density (e) divided by the temperature to the 
fourth power (T 4 ) for the cases of 2+1 and 2+1+1 flavors. 



2. The EOS at zero chemical potential 

Figure |] shows our results for the EOS with 2+1+1 flavors and compares it with the previously 
determined 2+1 flavor case [||]. The charm quark contribution grows with temperature, as expected 
and at the highest available T it contributes about 20% to the energy density. We conclude that in 
the cases where the charm quark is thermalized, its contribution to the EOS at temperatures higher 
than about 200 MeV, cannot be simply ignored. Our result is qualitatively similar to the previous 
work in Ref. |@|. 

3. The EOS at nonzero chemical potential 

Figure || presents our results for some of the pressure expansion coefficients which are directly 
related to the charm quark contribution at nonzero chemical potential. The first row shows the "un- 
mixed" coefficients and the second row - three of the "mixed" coefficients (they involve derivatives 
with respect to the light and heavy chemical potentials in addition to the charm one). The mixed 
coefficients are much noisier than the unmixed ones. We also find an unphysical behavior for the 
C002 coefficient, which around the transition region turns slightly negative. Possible causes could be 
the quenched approximation for the charm quark and the large heavy quark discretization effects. 
We intend to study this behavior in the near future. We also have results for the coefficients in the 
interaction measure expansion, but due to lack of space, we do not include them here. They will be 
published in a future work. 
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Figure 2: Some of the expansion coefficients for the pressure at nonzero chemical potential. 

The results for the contribution of the chemical potentials to the interaction measure and the 
energy density for the 2+1+1 flavor case are presented in Figs. |3| and ^ respectively. The results 
are for different values of jli/T , while \i^ c are tuned such that the heavy (strange) and the charm 
quark densities are zero within statistical errors. In these figures we also include our previous 2+1 
flavor results The comparison between these two cases shows a negligible contribution due to 
the charm quark addition at low temperatures. At high temperatures the statistical errors are too 
large for a definite conclusion about the charm quark effects. 



4. The isentropic EOS 

The heavy-ion collision experiments produce matter that expands isentropically. This implies 
that s/tib remains constant. For the data obtained at the AGS, SPS, and RHIC experiments, s/tib 
is approximately 30, 45 and 300, respectively. We have determined the approximate isentropic 
trajectories in the (;U/, jih, Ho T) space, by numerically solving the system 

S VI Yl 

—(jM,Hh,Pc)=C, -^(n h n h ,n c )=0, -^(fi h n h ,n c )=0, (4.1) 

riB r j r j 

with C = 30, 45 and 300. The results for the interaction measure, pressure and the energy density 
along these isentropic trajectories are shown in Figs. || and [7[ for both the 2+1+1 and the 
2+1 flavor cases |0J. The comparison indicates that the charm quark contribution is non-negligible, 
although it is due mostly to the contribution of the zeroth order coefficients in the Taylor expansions 
(i.e., the EOS calculated at zero chemical potential). 
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5. Conclusions 

The addition of the charm quark has a noticeable effects on the EOS at zero chemical potential 
at high temperatures. This implies that for the EOS of the early Universe we cannot ignore the 
charm quark contributions. It is possible that the rest of the flavors should be included too if high 
enough temperatures are studied. The effects of the charm quark on dl and de at nonzero chemical 
potential are small. Since the charm quark mass is large, many of the Taylor expansion coefficients 
which involve derivatives with respect to fl c /T do not reach the massless Stefan-Boltzmann values 
in the temperature range studied. We also have found a negative "dip" for C002 around the transition 
temperature — an unphysical effect whose origin we intend to study in the near future. And lastly, 
we found that the isentropic EOS is also affected by the addition of the charm quark (due mainly 
to its large effect on the EOS at zero chemical potential.) 
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Figure 4: The correction to the energy density due to the nonzero chemical potentials 
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Figure 5: The isentropic interaction measure normalized to the fourth power of temperature. 
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Figure 6: The isentropic pressure normalized to the fourth power of temperature. 
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Figure 7: The isentropic energy density normalized to the fourth power of temperature. 
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